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Question 1 [8 Marks] 


(a) A vector is input in MATLAB via the command: 
a.= [67452+123351% 


Write down the expected output from the following subsequent commands: 


(i) b = a.72 
(ii) c = a(a>=4) 
(iii) d = find(a<4) 
(iv) e = [a a] 

) 


(v) a(a>=1 & a<=3)=0 


(b) Consider the following function, func.m: 


function [f] = func(x,y) 
f=x72-2Q*x*yty~2; 


end 


What would be the associated output from the MATLAB command z=func (4,3)? 


(c) Assuming that x is a vector of real numbers, what does the following MATLAB function 
do? 


function [m,c] = func2(x) 
n = numel(x); 
m = x(1); 
c= 1; 
for ii = 2:n 
if xCii) <n 
m = x(ii); 
c = ii; 
end 
end 


end 
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Question 2 [8 Marks] 


(a) 


For a set of n numbers {21,%2,...,%n}, the geometric mean, m, is defined by: 


n \3 . 
m = (I>) = (Giioncstn) x 
i=1 


Write a MATLAB function to determine the geometric mean of all the elements of a 


matrix mat without using MATLAB’s built-in geomean command. 


Write a MATLAB script to evaluate y = x? — 9x and y’ = 3x? — 9 at 10000 equally 
spaced values of x over the interval —5 < x < 5, and then graph the functions as a 
red dashed line and a blue solid line, respectively. Add code to your script to label 
both the z- and y-axes of the graph, and add a legend to distinguish between the two 


Curves. 
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Question 3 [6 Marks] 


(a) Convert the following numbers from binary to decimal form, showing all working: 


(b) Convert the following numbers from decimal to binary, using the pen and paper algo- 


rithms demonstrated in lectures: 
(i) 24 

(ii) 0.75 

(iii) 24.75 


(c) Re-write the binary representation of 24.75 as a double precision floating point number 
(in binary) according to IEEE standard 754. 


(d) Consider a base 7 system for representing numbers that only uses the symbols 0, 1, 2, 


3, 4, 5 and 6. Convert the number represented by 10101, in binary to base 7. 


AMTH250 Exam Trimester 1 2019 


Question 4 [5 Marks] 


(a) 


The general Taylor series expansions for a function, f(a), expanded about a point 


X = Xp, and evaluated at the points x,4, and x,_, can be written as: 


h2 h3 hs. 
f(tau) = f(r) +hf' (tn) + aif (en) + af (en) + qf en) + O(h”) 


and 
/ h? " h? m hi (iv) 5 
f(@n-1) = flan) — hf" (en) + at (Zn) — 3r (fn) + a (Tn) + O(h?) 
respectively, where the points = ...,%n_1,%n,%n41,--- are equally spaced with sep- 


aration h = 24, —2;, and O(h°) indicates terms that are of the order of h® and above. 
Use the above Taylor series to obtain central difference approximations for f’(x,,) and 
f" (an) that have associated errors of O(h?). 


The following table lists outputs of some function f for some equally separated inputs 


x: 


f | 0.0000 | 0.7071 | 1.0000 | 0.7071 0.0000 
0.0000 | 0.7854 | 1.5708 | 2.3562 3.1416 | 


Use your central difference approximations from part (a) to estimate f’(2.3562) and 
f"(1.5708) based on the above table of values. 
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Question 5 |7 Marks] 
The following is an incomplete section of a MATLAB script to apply the bisection algorithm 


to find a zero of the function contained in a file func3.m within the interval [0, 10). 


clear 
x1=0; 


xr=10; 


f1=func3 (x1); 
while abs(xr-xl)>=... 4%Incomplete line 1 
xnew=... 4Incomplete line 2 


fnew=func3 (xnew) ; 


if ... %Incomplete line 3 

xl=xnew; 

fl=... %Incomplete line 4 
else 

xr=... %Incomplete line 5 
end 


end 


(a) Complete the indicated lines of code so that the zero is determined correct to 10 decimal 


places. 


(b) How many cycles of the while loop will it take to determine the zero to the desired 


level of accuracy? 
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Question 6 [10 Marks] 


2 


Consider the function g(a) =e” +1. 


(a) 


Sketch a simple graph of y = e~” +1 over the interval 0 < x < 2, along with a graph of 
y = x. Based on your graph, identify an approximate solution to x = g(x) that could 
be used as the initial input, xo, into a fixed point iteration algorithm to determine a 


more accurate solution. 


What is the general condition, in terms of g’(x), for fixed point iteration to converge to 
a solution of x = g(x)? Given that the root of x = e~®” +1 is approximately 1.2237, do 
you expect a fixed point iteration algorithm to converge to the solution of x = er 4 iF 


or not? 


Re-write the equation x = g(x) in the form f(x) = 0 (explicitly stating the exponential 
and constant terms). Then, write down the Newton-Raphson iteration formula for 


solving this specific problem. 


The following is an incomplete MATLAB script to determine the solution to the specific 


area using Newton-Raphson iteration. Complete the indicated lines 


equation 7 = e— 
of code so that the solution is determined to machine precision (using a relative error 


based stopping criterion), and with a maximum of 500 iterations. 


xold = 1.5; 
for ii = ... %Incomplete line 1 
xnew = ... %Incomplete line 2 
if abs(xnew-xold) < ... %Incomplete line 3 
break 
end 
. = xnew; %Incomplete line 4 


end 
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Question 7 [10 Marks] 
Consider the following data: 


aaa ae 2 3 4 
Graal Oe | 3.6 | 5 ip 
yi | 0} 0.84 | 0.6 | —0.44 | —0.96 


(a) Draw a sketch that illustrates the approximate positions of the given data points in 


the xy-plane. 


(b) What is the minimum degree you would expect for a polynomial that passes exactly 


through all five data points? 


(c) If you were to estimate y(2) using linear interpolation, which pair of points would be 


the most appropriate to use? 


(d) Write down the Lagrange polynomial for x73 = 5, that is, the polynomial L3(x) of 
smallest degree such that D3(%3) = 1 and Ls(a;) = 0 for all 7 € {0,1, 2, 4}. 


(ce) Assuming that L;(x),i = 0,...,4, are the Lagrange polynomials for the given x-values, 
write down the interpolating polynomial for the points given in the table as a linear 
combination of the polynomials L;(x). (There is no need to write out each of the 


polynomials, L;(x), explicitly for this question.) 


(f) Suppose, we want to construct a cubic spline f(a) through the given data points, that 
is, f is defined by f(x) = fi(x) for x € [2;, 241], 7 =0,1,...,3, where the functions f; 


are cubic polynomials which we can write in the form 
f(@) = ale —-2,) + b(2 = 2) +(e —2,) + di: 


For the specific data given in the table above, write down the conditions on the coef- 


ficients a;, b;, c,d; that come from the constraints 


(i) fo(v2) = ye (ii) fo(vs) = f3(73) 
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(g) What is an advantage of constructing the cubic spline using cubic polynomials of the 
form 


f(x) = a;(x — 2,)° + b;(x — 2)? + G(x — xj) + d; 


as opposed to cubic polynomials of the folowing form? 
fi(x) = a,x? + b,x" + Cx + d; 


(h) The MATLAB command, yy = spline(x,y,xx) returns a vector yy of values at the 
positions given in vector xx of a cubic spline interpolating the points (x(i),y(i)) for 
i=0, 1, ..., n. Write a short MATLAB script to plot the spline over the interval 
[(0, 10] and evaluate the spline at x = 4 for the data given in the table at the beginning 


of this question. 
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Question 8 [6 marks] 


From the values 


< 0 | 7/8 1/4 3n/8 | 1/2 
cos(x) | 1 | 0.9239 | 0.7071 | 0.3827 | 0 


a /2 
i: cos x dx 
0 


using the trapezoidal rule and Simpson’s rule, both with 3 and then 5 function evaluations. 


compute approximations for 


Determine the exact value of the integral and compare it to your approximations. 


Please remember: This examination question paper MUST BE HANDED IN. 


Failure to do so may result in the cancellation of all marks for this examination. 
Writing your name and number on the front will help us confirm that your paper 


has been returned. 
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